Generalized effective depletion potentials 
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We propose that the behavior of asymmetric binary fluid mixtures with a large class of attractive 
or repulsive interparticle interactions can be understood by mapping onto effective non-additive 
hard-sphere models. The latter are best analyzed in terms of their underlying depletion potentials 
which can be exactly scaled onto known additive ones. By tuning the non-additivity, a wide variety of 
attractive or repulsive generalized depletion potential shapes and associated phase behavior can be 
"engineered" , leading, for example, to two ways to stabilize colloidal suspensions by adding smaller 
particles. 
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Colloidal suspensions are complex mixtures of meso- 
scopic solute particles and smaller solvent particles. Ex- 
amples include a wide variety of organic or inorganic so- 
lutes, ranging from proteins to micelles to polymeric com- 
posites to ceramic materials etc., suspended in polar or 
non-polar solvents. Varying the interactions between the 
constituent particles results in a broad range of equilib- 
rium and non-equilibrium fluid behavior. This tunability 
has led to the wide spread industrial and biological appli- 
cations of colloidal suspensions . Some more recent de- 
velopments include the design of complex self-assembled 
materials such as photonic band-gap crystals by use of 
templates j^], and new experimental advances that al- 
low colloidal interactions to be directly measured with a 
greatly increased accuracy . 

In all the applications mentioned above, the design of 
a colloidal fluids with certain desired properties requires 
control over the inter particle interactions. These inter- 
actions are typically effective, that is to say they are 
a combination of direct interactions (such as Coulomb 
forces) with indirect interactions mediated through the 
solvent and the other solute particles |^-^ . One of the 
best known is the indirect depletion interaction, where 
one set of (typically smaller solute or solvent) particles 
induces an interaction between another set of particles. 
Depletion potentials were flrst calculated for mixtures of 
polymers and colloids and, with the advent of new 
experimental and theoretical techniques, they have been 
the subject of intense recent interest |p|-p^ 

Theoretical work has often focussed on the binary 
hard-sphere (HS) model, for which a depletion induced 
phase separation for size-ratios's q = ass/crbb < 0.2 was 
suggested [|| {aaa is the diameter of the big (6) or small 
(s) particles). A key advance was made by Dijkstra 
et al. JTof , who used an effective one-component deple- 
tion potential picture to show that the fluid-fluid phase- 
separation found with a two-component integral equation 
technique by Biben and Hansen [|j was metastable w.r.t. 
a fluid-solid phase-transition. More generally, their ap- 
proach added to the growing consensus that a carefully 



derived effective potential is a powerful tool for analyzing 
the behavior of an asymmetric binary mixture, at least 

for size-ratios q < 0.3 where many-body interactions are 
not thought to be important (see e.g. for some re- 

cent reviews). The key step in all these approaches is in- 
tegrating out the smaller component of a binary mixture 
to leave a new one-component system with an effective 
interaction between the big particles. 

Most theories of depletion have considered only hard- 
core interactions leading to purely entropic depletion 
potentials. Their range varies with ass while increas- 
ing the small particle density ps or packing fraction 
77s = TTpsCTss/Q increases the depth of the (always) at- 
tractive well at contact, and possibly adds enhanced os- 
cillations at larger separations r [@j|,|l^ . 

There have been a number of recent attempts to go 
beyond purely entropic depletion by including extra in- 
teractions between the particles of a binary HS mixture 
Ip^-p^. Of course many different kinds of extra inter- 
actions can be added, leading to a seemingly enormous 
increase in complexity. However, in this letter we pro- 
pose that the effect of a wide variety of these extra inter- 
actions on depletion potentials can be understood by a 
simple mapping onto a non-additive HS mixture model, 
for which the depletion potentials can be calculated by a 
second exact mapping or scaling onto those of an additive 
system. 

Since the phase behavior of many binary fluids can be 
well understood on the basis of these depletion potentials 
iQ-^ , this implies that our (double) mapping can be used 
to analyze a wide variety of interacting asymmetric bi- 
nary mixtures. These ideas can also be turned around, 
leading to the possibility of explicitly engineering a wide 
variety of generalized depletion potential shapes, includ- 
ing potentials that are repulsive at contact, by not only 
varying the usual parameters ps and q, but also by tuning 
the interparticle interactions to vary the non-additivity. 

Non-additive binary HS models are deflned by specify- 
ing the cross-diameter [pOl 
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(Jbs = ^{(Jss+<7bb) (l + A). (1) 

When A = 0, the model foUows the Lorentz mixing rule, 
and is traditionally called additive (not to be confused 
with pairwise additivity of potentials), that is the cross- 
diameter is simply the sum of the two radii, exactly 
what one would expect on purely geometric grounds. 
If A > or A < the system shows positive or neg- 
ative non- additivity respectively. As shown in Fig. |^, 
each big particle excludes a volume Vb = from 
the centers of the smaller particles. When the deple- 
tion layers of the two big particles (width defined as 
h = abs-crbbfi = i(crss + A(CTss+crb6))) begin to overlap, 
then the small particles can gain free volume ua, leading 
to a depletion interaction. 



A>0 A<0 




FIG. 1. Centers of the small-particles of diameter asa 
are excluded from the shaded depletion layer of width h 
around each big particle of diameter abb- If /i < fcTss 
then A < (negative non-additivity) . If /i > ^o-ss, then 
A > (positive non-additivity) 

To calculate these potentials we first note that the 
depletion potential pVeffir) depends only on the big- 
small and small-small interactions {pVbs{r) and pVss{r) 
respectively), but not on any direct big-big interaction 
/3Vbb(r), which can simply be added to the depletion po- 
tential [ppr^. For non-additive systems at fixed ps this 
means that the depletion potential is determined by <7bs 
and <7ss, and is equivalent to an additive one with the 
same parameters! Only the cutoff due to a^b is differ- 
ent. For example, if one has an expression for the addi- 
tive potential in terms of pb, q and the scaled distance 
p = r/abb, then the a potential with A 7^ is given by: 
f3Veff{ps,q,A,p) = pVeff{ps,q\^ = 0,(<zVg)p) where 
q' = ass/{'2crbs — Gss/'^)- Details of this (perhaps sur- 
prising) exact mapping by scaling to additive potentials 
will be given elsewhere |2|] . The depletion potentials for 
additive systems can be calculated to quantitative ac- 
curacy by a density functional theory (DFT) technique 
, and we checked that the scaling procedure above ex- 
actly reproduces recent DFT calculations of non-additive 
depletion potentials JTot . 

For a fixed number density ps, non-additivity can be 
introduced in two ways: 

Case (A) Fix the depletion layer width h, (or equiv- 
alently abs), and vary A by changing the small parti- 
cle diameter Gss- The effect on depletion pair-potentials 
j3Veff{r) is shown in Fig. ^. For increasing positive non- 
additivity the correlation induced maximum decreases 



and the potential tends towards the (ideal) Asakura- 
Oosawa (AO) limit; the contact value remains rela- 
tively constant, as was found earlier In contrast, 
for increasing negative non-additivity the contact value 
increases markedly, leading to the possibility of strongly 
repulsive interactions. A naive application of the simplest 
depletion picture where (3Veff(r = abb) = — H^wa would 
give the opposite effect, since decreasing A by increas- 
ing ass increases the packing fraction rjs — irpsass/Q and 
therefore the small particle osmotic pressure Us, while 
keeping va unchanged, seemingly leading to a more at- 
tractive effective potential. However, a more careful anal- 
ysis reveals that increasing 77^ leads to well-developed 
solvation shells around a single big particle. When two 
big particles approach, the overlap of the solvation shells 
leads to the repulsive interactions, as well as larger os- 
cillations in the pair potentials. We found that the am- 
plitude of this repulsion becomes larger for smaller size- 
ratios, with values possible of many times fcsT. 

Case (B) Fix the small-particle hard-core diameter 
ass , and vary A by changing h (or equivalently abs ) ■ The 
dominant effect on depletion pair-potentials is to shift 
them along r as shown in Fig. |3[ In this case both pos- 
itive and negative non-additivity change the well-depth 
at contact significantly because changing h changes the 
amount of volume doubly excluded when two big parti- 
cles approach. This can be illustrated at the simple AO 
level where the potential at contact is given by 

pVAoir - abb) - -Ps^ (^abb{2h)^ + ^(2^') ■ (2) 

On the other hand, the correlation induced maximum re- 
mains roughly the same since rjs is constant, leading to 
similar solvation layers of the small particles around a 
big particle. 
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FIG. 2. Depletion potentials when parameters are 
changed according to case (A), i.e. varying Gss but keep- 
ing h constant. Here q — 0.2, and A-Kpgh^ /i = O.I. The 
potentials are calculated by an exact scaling from known 
additive DFT results 

With the insight gained from analyzing non-additive 
depletion potentials, we next pose the question: what 
happens to /3T4//(r) for a given binary HS mixture 
when more general attractive or repulsive interactions 
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f3Vss{r) or Vbs{r) are added? Inspired by some well- 
established ideas from the theory of simple liquids [ pO| , 
we map onto effective HS diameters as follows: ctq,^ = 
+ /(6xp[— /3yQ/3(r)] — l)dr, a procedure similar to 
the well known Barker-Henderson approach [^. Here 
cr°^ denotes the original effective HS diameters without 
the extra interaction. In this way the additional interac- 
tions can be mapped onto an effective non-additive HS 
model as follows: 

,0 . „ ^ ^0 



(i) repulsive l3Vss{r) 

(ii) attractive pVss{r) : 

(iii) repulsive f3Vbs{r) 

(iv) attractive (3Vbs{r) 
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FIG. 3. Depletion potentials when parameters are 
changed according to Case (B), i.e. varying h but keep- 
ing ass constant. Here q = 0.2, rjs — 0.2 and line-styles 
denote the same As as in Fig. ^. Inset: Depletion po- 
tentials when l3Vbs{r) — eexp[—K,{r — abs)]/r is added to 
a binary HS mixture are compared to the non-additive 
HS case with A calculated by our simple mapping. Here 
rjs = 0.116, 5 = 0.2, KCTss = 4 is fixed and f3e is varied. 
The included simulation data for jSe = —0.32 |^ helps 
confirm the accuracy of our new direct DFT approach. 

The depletion potentials for (i) (A < 0) and (ii) 
(A > 0) change according to case (A), which is depicted 
in Fig. 1^, while the depletion potentials for (iii) (A > 0) 
and (iv) (A < 0) change according to case (B), as de- 
picted in Fig. ^. This picture agrees qualitatively with 
calculations of other authors of depletion potentials for 
non-HS systems. Examples of the pathways above in- 
clude: (i) Fig. 4 of |18|; (ii) Fig. 8 of pf. F ig. 4 of §; 
(iii) Fig. 5 of |T4|; Fig. 3 of |l^; (See also||3| for a recent 



experimental and [p4| for a recent theoretical mapping 
of this type of interaction onto an AO potential), (iv) 
Fig. 4 of |l|]. Fig. 4 of Hi, Fig. 7 of |l7|, and Fig. 4 of [|. 
Keep in mind, however, that these calculations were done 
with a number of approximate techniques which may not 
always give quantitatively reliable results, especially for 
contact values Q. Even so, our mapping scheme to non- 
additivity clearly qualitatively rationalizes the dominant 
changes in depletion potentials caused by changing a va- 
riety of inter-particle interactions. 

It would be interesting to make this qualitative cor- 
respondence more quantitative. We were able to extend 



the quantitatively reliable DFT method described in [T^ ] 
to systems with an arbitrary Vbs (r) potential (details will 
be published elsewhere). In the inset of Fig ^ we compare 
these to non-additive HS potentials with A determined 
by our aforementioned mapping procedure. This gives 
a reasonably good representation of the well depth, but 
does slightly less well for the repulsive barrier. With 
our very simple mapping procedure, we expect that the 
quantitative agreement will deteriorate for very strongly 
attractive Vbs{r) or Vss{r), but the qualitative picture 
should remain the same. 

From the above it is clear that non-additivity has a 
profound effect on depletion potentials, implying that 
this should also be reflected in phase-behavior. For more 
symmetric mixtures the effect of interparticle interactions 
on phase-stability has traditionally been understood in 
terms of conformal-solution theory ||20|| . For asymmetric 
mixtures, our depletion potentials help generalize these 
ideas. For example, Vliegenthart and Lekkerkerker ]2^ ] 
have recently shown that the fluid-fluid critical point of 
many one-component fluids occurs when the reduced sec- 
ond virial-coefficient 82/62^ w —1.5. We checked that 
this works well for the depletion potential simulations of 
Dijkstra et al. [ poHTl| , suggesting that this surprisingly 
accurate criterion can also be used to predict the effect 
of non-additivity on fluid-fluid phase-separation. In Fig ^ 
we plot the effect of the rjs on the second virial coefficients 
calculated from depletion potentials. Firstly, for the ad- 
ditive case we find that 62/82^ < —1.5 only for size- 
ratios q < 0.11 (For q = 0.1 there might be an upper crit- 
ical point!). But even without attributing quantitative 
accuracy to the Vliegenthart Lekkerkerker criterion, the 
upper limit of q that allows fluid-fluid phase-separation is 
clearly bounded by 0.1 < q < 0.2, since 62/62^ remains 
positive for any 77^ if g > 0.2, while it goes well below 
— 1.5 for q < 0.1. This finding helps in understanding 
earlier results obtained with (approximate) 2-component 
integral equation studies as well as some direct sim- 
ulations |]l^ , lending support to our argument that the 
underlying depletion potentials from which the B2 are de- 
rived are a key to understanding the full phase-behavior 
of the asymmetric two-component systems. 

Next we turn to the effect of non-additivity on the 
fluid-fiuid phase-separation. Fig. || shows that for g = 0.1 
a very small non-additivity, of the order of a 5% change 
in (Jss or a 0.5% change in abs, is enough to dramatically 
change the behavior of B2/B2^ . For other size-ratios 
we find similar effects. For example if A = (7/2O wc find 
that (metastable) fluid-fluid phase-separation can occur 
for size-ratios up to 5 = 0.4, while if A = —q/20, it 
will only occur for size-ratios q < 0.05. Clearly, even 
a very small negative non-additivity strongly suppresses 
phase-separation, while positive non-additivity strongly 
enhances it. This is consistent with and helps rationalize 
some earlier 2-component studies [p6|-p8[. 

Binary mixtures may also undergo fluid-solid phase- 
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separation which, for example, is the thermodynamically 
stable phase-transition in additive HS mixtures | [To[ | . Re- 
cently, one of us |^,^ has shown that for short range 
potentials the fluid-solid transition shifts to low values 
of rjb = 7r/9{,ajfj/6 when the potential well-depth at con- 
tact is near (jVeff{r ~ abb) ~ —2.5; this effect is largely 
independent of other details such as the range or oscil- 
lations of the potential. This suggests that introducing 
any non-additivity according to case (B) will strongly 
affect the fluid-solid behavior. Similarly adding negative 
non-additivity according to case (A) will suppress fluid- 
solid phase-separation, but positive non-additivity will 
not change the fluid-solid phase-boundaries much, which 
is confirmed by comparing the additive (A = 0) HS to 
pure AO (q — A) simulations of Dijkstra et al. |lo| , [rT|] . 
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FIG. 4. Reduced second virial coefficient of the big 
particles B2/B^^, plotted v.s. the packing fraction of 
the small particles rjs . Fluid-fluid phase-separation is ex- 
pected when B2/B2^ < -1-5 (25) (horizontal line). 

The stabilization of colloidal suspensions is critical 
to many industrial and experimental applications [Q. 
The arguments above for both fluid-fluid and fluid-solid 
phase-separation suggest that the addition of smaller par- 
ticles may provide such a stabilization mechanism against 
demixing for colloidal suspensions as long as A < 0. 
This can be achieved by pathway (i) , adding a repulsive 
Vss{r), or by pathway (iv), adding an attractive Vbs(r). 

In conclusion then, we combined a new approximate 
and a new exact mapping to show that the non-additive 
HS mixture model provides an intuitive and general or- 
ganizing framework within which to understand the ef- 
fective depletion potentials induced by a large class of 
interactions Vbs{r) or Vssir). These generalized deple- 
tion potentials can be crafted into many different shapes, 
and provide access to the phase-behavior of interacting 
asymmetric binary mixtures. Clearly much more can be 
done by both theories and experiments to exploit the 
flexibility of these potentials and to "engineer" desired 
phase-behavior in colloidal suspensions. We hope this 
letter has shown some promising new directions in which 
to embark. 
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